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Abstract 

We discuss the current understanding of nuclear scattering at very high energies. We point 
out several serious inconsistencies in nowadays models, which provide big problems for any 
interpretation of data at high energy nuclear collisions. 

We outline how to develop a fully self-consistent formalism, which in addition uses all the 
knowledge available from studying electron-positron annihilation and deep inelastic scattering, 
providing a solid basis for further developments concerning secondary interactions. 



1 Introduction 

We will report on efforts to construct realistic and reliable models for nuclear collisions at very 
high energies, say above 100 GeV cms energy per nucleon. There are several reasons why such 
models are needed, related to large scale experimental programs presently and in the future. 

In astrophysics, one of the open problems is the origin of cosmic rays at very high energies. 
Whereas at low energies direct measurements are possible using balloons or satellites, at high 
energies - due to the small flux - only indirect measurements are possible, counting the secondary 
particles created in a so-called extended air-shower. Here, theoretical models for hadron-nucleus 
and nucleus-nucleus collisions play a decisive role in order to provide some information on the 
origin of the primary particle. 

In nuclear physics, there are ongoing efforts to produce and to investigate the so-called quark- 
gluon plasma, a new state of matter of quarks and gluons. The SPS program at CERN has revealed 
many interesting results, but a final conclusion will never be made without a profound theoretical 
understanding, based on realistic microscopic models. 

The outline of this paper is as follows: we will first discuss the present status concerning 
models for the initial stage of a nucleus-nucleus coUision. We find that even the most advanced 
approaches show severe theoretical inconsistencies, which make any interpretation of experimental 
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data questionable. It is in fact a well known problem which has been brought up a decade ago, but 
no solution has been proposed so far. We will discuss how to solve the above-mentioned problem, 
introducing a self-consistent formalism for nuclear scattering at very high energies. 

2 Present Status 

Many popular models 0, ||, [sj are based on the so-called Gribov-Regge theory (GRT) |^, ^. This 
is an effective field theory, which allows multiple interactions to happen "in parallel", with the phe- 
nomenological object called "Pomeron" representing an elementary interaction. Using the general 
rules of field theory, one may express cross sections in terms of a couple of parameters character- 
izing the Pomeron. Interference terms are crucial, they assure the unitarity of the theory. Here 
one observes an inconsistency: the fact that energy needs to be shared between many Pomerons in 
case of multiple scattering is well taken into account when calculating particle production (in par- 
ticular in Monte Carlo applications) , but energy conservation is not taken care of in cross section 
calculations. This is a serious problem and makes the whole approach inconsistent. 

Provided factorization works for nuclear collisions, one may employ the parton model |^, 0|, 
which allows to calculate inclusive cross sections as a convolution of an elementary cross section 
with parton distribution functions, with these distribution functions taken from deep inelastic 
scattering. In order to get exclusive parton level cross sections, some additional assumptions are 
needed, which follow quite closely the Gribov-Regge approach, encountering the same difficulties. 

Before presenting new theoretical ideas, we want to discuss the open problems in the parton 
model approach and in Gribov-Regge theory somewhat more in detail. 

Gribov-Regge Theory 

Gribov-Regge theory is by construction a multiple scattering theory. The elementary interactions 
are realized by complex objects called "Pomerons", who's precise nature is not known, and which 
are therefore simply parameterized: the elastic amplitude T corresponding to a single Pomeron 
exchange is given as 

T{s,t) - js°«+"'* 

with a couple of parameters to be determined by experiment. Even in hadron-hadron scattering, 
several of these Pomerons are exchanged in parallel. Using general rules of field theory (cutting 
rules), one obtains an expression for the inelastic cross section, 

'^it^(«)= / rf'Ml-cxp(-G(s,6))} = ^a:^^''^(.), (1) 

with 

^t'^^s) = f d25Mfi^exp(~G(.,6)), (2) 
J m\ 

where the so-called eikonal G{s,b) (proportional to the Fourier transform of T{s,t)) represents 
one elementary interaction. The cross sections al^''-^{s) for m inelastic collisions are referred to as 
topological cross sections. One can generalize to nucleus-nucleus collisions, where corresponding 
formulas for cross sections may be derived. 

In order to calculate exclusive particle production, one needs to know how to share the energy 
between the individual elementary interactions in case of multiple scattering. We do not want to 
discuss the different recipes used to do the energy sharing (in particular in Monte Carlo applica- 
tions). The point is, whatever procedure is used, this is not taken into account in the calculation of 
cross sections discussed above. So, actually, one is using two different models for cross section cal- 
culations and for treating particle production. Taking energy conservation into account in exactly 
the same way will modify the cross section results considerably. 
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This problem has first been discussed in The authors claim that following from the non- 

planar structure of the corresponding diagrams, conserving energy and momentum in a consistent 
way is crucial, and therefore the incident energy has to be shared between the diflerent elementary 
interactions, both real and virtual ones. 

Another very unpleasant and unsatisfactory feature of most "recipes" for particle production is 
the fact, that the second Pomeron and the subsequent ones are treated differently than the first 
one, although in the above-mentioned formula for the cross sections all Pomerons are considered 
to be identical. 

The Parton Model 

The standard parton model approach to hadron-hadron or also nucleus-nucleus scattering amounts 
to presenting the partons of projectile and target by momentum distribution functions, fhi and 
fh2, and calculating inclusive cross sections for the production of parton jets with the squared 
transverse momentum larger than some cutoff Qq as 

- E/ '^pIJ d^^J dx-fl^{x\pl)fi^ix-,pl)^ix+x-s)e{pl - Ql) , 

where daij/dp\ is the elementary parton-parton cross section and i,i represent parton flavors. 

This simple factorization formula is the result of cancelations of complicated diagrams (AGK 
cancelations) and hides therefore the complicated multiple scattering structure of the reaction. 
The most obvious manifestation of such a structure is the fact that at high energies (^/s ^ 10 
GeV) the inclusive cross section in proton-(anti-)proton scattering exceeds the total one, so the 
average number Nf^^ of elementary interactions must be greater than one: 

AThih2 hih2 /hih2 ^ -i 

^^int — "incl /"tot ^ ^ ■ 

The usual solution is the so-called eikonalization, which amounts to re-introducing multiple scat- 
tering, based on the above formula for the inclusive cross section: 

^it^(^) = / rf'fe {l-exp(-A(6)a[:^t^(.))} ^J^a:;,^''^^), (3) 

with 

/■ f^(^)^i',\cr(5))" / X 

at'-^i.) = j ^T^^^P [-Mh)at^,^{s)) (4) 

representing the cross section for n scatterings; A{b) being the proton-proton overlap function 
(the convolution of two proton profiles) . In this way the multiple scattering is "recovered". This 
makes the approach formally equivalent to Gribov-Regge theory, encountering therefore the same 
problems: energy conservation is not at all taken care of in the above formulas for cross section 
calculations. 

3 A Solution: Parton-based Gribov-Regge Theory 

In this paper, we present a new approach for hadronic interactions and for the initial stage of 
nuclear collisions, which is able to solve the above-mentioned problems. We provide a rigorous 
treatment of the multiple scattering aspect, such that questions as energy conservation are clearly 
determined by the rules of field theory, both for cross section and particle production calculations. 
In both (!) cases, energy is properly shared between the different interactions happening in parallel, 
see fig. 1^ for proton-proton and fig. ^for proton-nucleus coUisions (generalization to nucleus- nucleus 
is obvious). This is the most important and new aspect of our approach, which we consider to be 
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Figure 1: Graphical representation of a contribution to the elastic amplitude of proton-proton scattering. Here, 
energy conservation is taken into account: the energy of the incoming protons is shared among several "constituents" 
(shown by splitting the nucleon lines into several constituent lines), and so each Pomeron disposes of only a fraction 
of the total energy, such that the total energy is conserved. 




Figure 2: Graphical representation of a contribution to the elastic amplitude of proton-nucleus scattering, or more 
precisely a proton interacting with (for simplicity) two target nucleons, taking into account energy conservation. 
Here, the energy of the incoming proton is shared between all the constituents, which now provide the energy for 
interacting with two target nucleons. 

a first necessary step to construct a consistent model for high energy nuclear scattering. 

The elementary interactions, shown as the thick lines in the above figures, are in fact a sum 
of a soft, a hard, and a semi-hard contribution, providing a consistent treatment of soft and hard 
scattering. To some extend, our approach provides a link between the Gribov-Regge approach and 
the parton model, we call it "Parton-based Gribov-Regge Theory". 

4 Parton-Parton Scattering 

Let us first investigate parton-parton scattering, before constructing a multiple scattering theory 
for hadronic and nuclear scattering. 

We distinguish three types of elementary parton-parton scatterings, referred to as "soft", "hard" 
and "semi-hard", which we are going to discuss briefiy in the following. The detailed derivations 
can be found in |[ic|| . 

The Soft Contribution 

Let us first consider the pure non-perturbative contribution, where all virtual partons appearing in 
the internal structure of the diagram have restricted virtualities <Q^, where Qq ~ 1 GeV^ is a 
reasonable cutoff for perturbative QCD being applicable. Such soft non-perturbative dynamics is 
known to dominate hadron-hadron interactions at not too high energies. Lacking methods to cal- 
culate this contribution from first principles, it is simply parameterized and graphically represented 
as a 'blob', see fig. ^. It is traditionally assumed to correspond to multi-peripheral production of 
partons (and final hadrons) [Q and is described by the phenomenological soft Pomeron exchange 
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Figure 3: The soft contribution. 



amplitude Tsoh{s,t) |^. The corresponding profile function is twice the imaginary part of the 
Fourier transform Tsoft of Tgoft, divided by the initial parton fiux 2s, 




with 



^soft(^) = "-Rpart + "soft InZ 



where s is the usual Mandelstam variable for parton-parton scattering. The parameters asoft(O), 
a soft are the intercept and the slope of the Pomeron trajectory, 7part and i?part ^^'6 the vertex value 
and the slope for the Pomeron-parton coupling, and so — 1 GeV^ is the characteristic hadronic 
mass scale. The external legs of the diagram of fig. ||are "partonic constituents", which are assumed 
to be quark-anti-quark pairs. 

The Hard Contribution 

Let us now consider the other extreme, when all the processes are perturbative, i.e. all internal in- 
termediate partons are characterized by large virtualities > Qq. In that case, the corresponding 
hard parton-parton scattering amplitude can be calculated using the perturbative QCD techniques 
[T^ , and the intermediate states contributing to the absorptive part of the amplitude can 
be defined in the parton basis. In the leading logarithmic approximation of QCD, summing up 
terms where each (small) running QCD coupling constant as{Q^) appears together with a large 
logarithm \n{Q'^ / Xq^^)) (with Xqcd being the infrared QCD scale), and making use of the factor- 
ization hypothesis, one obtains the contribution of the corresponding cut diagram for t = = 
as the cut parton ladder cross section crj^ard('^' Qo) B ^® shown in fig. ^, where all horizontal rungs 




Bom 



QCD 



Figure 4: Ttie iiard (or val-val) contribution. 
^ Strictly speaking, one obtains ttie ladder representation for the process only using axial gauge. 
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are the final (on-shell) partons and the virtuahties of the virtual t-channel partons increase from 
the ends of the ladder towards the largest momentum transfer parton-parton process (indicated 
symbolically by the 'blob' in the middle of the ladder): 

<ard(S'<5o) 



X 

Here da'^l^^/ dp\ is the diflerential 2^2 parton scattering cross section, p\ is the parton transverse 
momentum in the hard process, m, I and are correspondingly the types and the shares of the 
Ught cone momenta of the partons participating in the hard process, and Mp is the factorization 
scale for the process (we use Mp = p\/'^). The 'evolution function' i?Q™j-,((5Q, M|., z) represents 
the evolution of a parton cascade from the scale Qo to Mp, i.e. it gives the number density of 
partons of type m with the momentum share z at the virtuality scale M|,, resulted from the 
evolution of the initial parton j' taken at the virtuality scale Qq- The evolution function satisfies 
the usual DGLAP equation [Q with the initial condition i?Q™j-,((5§, Qq, z) — d-j^ (5(1 — z). The 
factor K ~ 1.5 takes effectively into account higher order QCD corrections. 

In the following, we shall need to know the contribution of the uncut parton ladder tI^^^{s, t, Qq) 
with some momentum transfer q along the ladder (with t = q^). The behavior of the corresponding 
amplitudes was studied in ||l5|| in the leading logarithmic (1 /a; ) approximation of QCD. The precise 
form of the corresponding amplitude is not important for our appHcation; we just use some of the 
results of jl^, namely that one can neglect the real part of this amplitude and that it is nearly 
independent on i, i.e. that the slope of the hard interaction i?hard neghgible small, i.e. compared 
to the soft Pomeron slope one has iJ^ard — 0- So we parameterize 7haid('^' *3o) i'^ the region of 
small t as 

^hard(«' Ql) = f^hard(«' QD ^XP (i?hard (6) 

The corresponding profile function is twice the imaginary part of the Fourier transform Thard 
of Thard, divided by the initial parton flux 2s, 

which gives 



^2Im7]^,^,,(s,t = 0,Q^) 

/J m/ 

nil 

^qcd(Qo, Ml, 4) 4'cd(Qo, Ml, Xs)9{Ml - Ql) , 



^^^hard \ ^-^^hard / 

In fact, the above considerations are only correct for valence quarks, as discussed in detail in 
the next section. Therefore, we also talk about "valence-valence" contribution and we use i^vai-vai 
instead of Uhard: 

so these are two names for one and the same object. 
The Semi-hard Contribution 



The discussion of the preceding section is not vaHd in case of sea quarks and gluons, since here 
the momentum share xi of the "flrst" parton is typically very small, leading to an object with 
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Figure 5: The semi-hard "sea-sea" contribution: parton ladder plus "soft ends". 

a large mass of the order Qq/xi between the parton and the proton jl^. Microscopically, such 
'slow' partons with a;i <C 1 appear as a result of a long non-perturbative parton cascade, where 
each individual parton branching is characterized by a small momentum transfer squared < Qq 
1^, H^. When calculating proton structure functions or high-pt jet production cross sections this 
non-perturbative contribution is usually included in parameterized initial parton momentum distri- 
butions at — Qq- However, the description of inelastic hadronic interactions requires to treat it 
explicitly in order to account for secondary particles produced during such non-perturbative parton 
pre-evolution, and to describe correctly energy-momentum sharing between multiple elementary 
scatterings. As the underlying dynamics appears to be identical to the one of soft parton-parton 
scattering considered above, we treat this soft pre-evolution as the usual soft Pomeron emission, 
as discussed in detail in | [lO| . 

So for sea quarks and gluons, we consider so-called semi-hard interactions between parton 
constituents of initial hadrons, represented by a parton ladder with "soft ends", see fig. |[ As in 
the case of soft scattering, the external legs are quark-anti-quark pairs, connected to soft Pomerons. 
The outer partons of the ladder are on both sides sea quarks or gluons (therefore the index "sea- 
sea"). The central part is exactly the hard scattering considered in the preceding section. As 
discussed in length in the mathematical expression for the corresponding amplitude is given 
as 

jk 

with being the momentum fraction of the external leg-partons of the parton ladder relative to 
the momenta of the initial (constituent) partons. The indices j and k refer to the flavor of these 
external ladder partons. The amplitudes T^^f^ are the soft Pomeron amplitudes discussed earlier, 
but with modified couplings, since the Pomerons are now connected to a parton ladder on one side. 
The arguments sq/z^ are the squared masses of the two soft Pomerons, z~^z~s is the squared mass 
of the hard piece. 

Performing as usual the Fourier transform to the impact parameter representation and dividing 
by 2s, we obtain the profile function 

-^sca— sea (^i ^) 77^ 2Im T^gca— sca('^, ^) , 
ZS 

which may be written as 
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("2") ' h 

with the soft Pomeron slope A^^j^ and the cross section <J^^^.^ being defined earUer. The functions 
-^soft i^^) representing the "soft ends" are defined as 

We neglected the small hard scattering slope -R^ard compared to the Pomeron slope A^^lj. We call 
E'soft also the " soft evolution", to indicate that we consider this as simply a continuation of the 
QCD evolution, however, in a region where perturbative techniques do not apply any more. As 
discussed in |[lo), E^^^^ (z) has the meaning of the momentum distribution of parton j in the soft 
Pomeron. 

Consistency requires to also consider the mixed semi-hard contributions with a valence quark 
on one side and a non-valence participant (quark-anti-quark pair) on the other one, see fig. §. We 




Figure 6: Two "mixed" contributions. 

have 

^T^4l-sea(^) = E I-^s'^.i^, l') i^-S, l\Ql) 

and 

DL-soS^M = J2 t'^'^Esoui^-X'U^-S^Ql) (9) 

1 f 9 

^ — m — m 

4^aW (1/Z-) V 4A« (1/Z-) 

where j is the flavor of the valence quark at the upper end of the ladder and k is the type of the 
parton on the lower ladder end. Again, we neglected the hard scattering slope i^h^rd compared to 
the soft Pomeron slope. A contribution fgoa-vai corresponding to a valence quark participant 

from the target hadron, is given by the same expression, 

^sca-val (S' ^) = ^val-soa ^) ' 

since eq. (||) stays unchanged under replacement z+ and only depends on the total cm. 

energy squared s for the parton-parton system. 



5 Hadron-Hadron Scattering 

To treat hadron-hadron scattering we use parton momentum Fock state expansion of hadron eigen- 
states 1^ 

oc „i A: / ^' \ 

1^) = E / n^^' (^1' ■ • ■ ^fc) '5 1 - E^^ • • • io) ' 
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where fk{xi, . . . Xk) is the probabiHty amphtude for the hadron h to consist of k constituent 
partons with the Hght cone momentum fractions xi, . . . ,Xk and a+ [x) is the creation operator for 
a parton with the fraction x. A general scattering process is described as a superposition of a 
number of pair-Hke scatterings between parton constituents of the projectile and target hadrons. 
Then hadron-hadron scattering amplitude is obtained as a convolution of individual parton-parton 
scattering amplitudes considered in the previous section and "inclusive" momentum distributions 
~\^h^^ (xi, . . - Xn) of n "participating" parton constituents involved in the scattering process(n > 1), 
with 

-f^"^(xi,...a;„) = I]^ ^,(fel^), l[dxi \fk{xi,...Xktsh-J2^j 

We assume that -^/^"(/ij) (^i' ■ ■ •^") '-^^ represented in a factorized form as a product of the 
contributions F^^^^{xi), depending on the momentum shares xi of the "participating" or "active" 
parton constituents, and on the function -Frcmn(l ^ Sj^i ^i) i representing the contribution of all 

spectators 

nucleon 

participants: 1 2 ... k 




Figure 7: Nucleon Fock state. 

"spectator" partons, sharing the remaining share 1 — Ylj of the initial Hght cone momentum (see 
fig- I): 

n / ^ \ 

Fi") {x^,... Xn) = n^p",,, {xi) h'J2^A (10) 

The participating parton constituents are assumed to be quark-anti-quark pairs (not necessarily of 
identical flavors), such that the baryon numbers of the projectile and of the target are conserved. 
Then, as discussed in detail in ||l^, the hadron-hadron amplitude may be written as 

tn,hAs,t)^87r^sY,- / Y[dx+dx; H d'qi^iT^^''-{x+,x-,s,-ql) 
11=1 1=1 i=i L 

^^rcU|^l - j ^^rt.„|^l - E^7 j ^^'^ (E-^^^ - ' (11) 

where the partonic amplitudes are deflned as 

rphih'^ rr\hih2 i rr\hih2 i rr\hih2 j, rr\hih2 i rphih2 

IP soft ~'~ sea— sea ' val — val ' val — sea ' sea— val' 

with the individual contributions representing the "elementary partonic interactions plus external 
legs". The soft or semi-hard sea-sea contributions are given as 

^stft'/soa-seal^''''^"''*'"?!) = T'soft/sea-scals, -?!) ^part (^^'^ ) -^part (2^" ) 

y. e^v{-[Rl^+Ri;\ql), (12) 
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the hard contribution is 

Jo Xvi Jo Xv . J, 

X F^l;iixt,x+ - xt)F!;i^{x-,x- - X-) exp(- [rI^ + i^^J ^fj ^ 
the mixed semi-hard "val-sea" contribution is given as 

Ttl'-U^^^-,^,^) - f dxt"-Y^n^,_^^^{xtx-s,q\Ql) 

Jo "^"^ j 

X F^l;({xt,x+ - xt)F^i,,{x-) exp(- [rI + R^] qfj , 
and the contribution "sea-val" is finally obtained from "val-sea" by exchanging variables, 

Here, we allow formally any number of valence type interactions (based on the fact that multiple 
valence type processes give negligible contribution). In the valence contributions, we have convo- 
lutions of hard parton scattering amplitudes T^^rd ^^"^ valence quark distributions -Fpa^t ^^^^ the 
valence quark momentum share x^ rather than a simple product, since only the valence quarks are 
involved in the interactions, with the anti-quarks staying idle (the external legs carrying momenta 
x~^ and x" are always quark-anti-quark pairs). 
The profile function 7 is as usual defined as 

lhih2{s,b) = ^2Imf ^1/12(5,6), 

which may be evaluated using the AGK cutting rules with the result (assuming imaginary ampli- 
tudes) 



CO ^ ^1 m m 

Hi. ft) 

m=l •''^ tj.= l /j=l 

00 „i / / 

Eli / l[dS:tdi-,l[-Gl^^^{xt,x-,,s,b) 

V A / V A / 

with xP™^/*'^''^ = 1 ~ S^/i being the momentum fraction of the projectile/target remnant, and 
with a partonic profile function G given as 

Gl^''^{x+,x^,.s,b) = j-2lmf^^^^ixt,x^,.s,bl (14) 



see fig. 1^. This is a very important result, allowing to express the total profile function 7/11/12 via 

'ITP 



the elementary profile functions G\^^^^ 



6 Nucleus-Nucleus Scattering 

We generalize the discussion of the last section in order to treat nucleus-nucleus scattering. In 
the Glauber-Gribov approach ||l^, ||, the nucleus-nucleus scattering amplitude is defined by the 
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hih-2 



Figure 8: The hadronic profile function 7 expressed in terms of partonic profile functions G = 



sum of contributions of diagrams, corresponding to multiple scattering processes between parton 
constituents of projectile and target nucleons. Nuclear form factors are supposed to be defined by 
the nuclear ground state wave functions. Assuming the nucleons to be uncorrelated, one can make 
the Fourier transform to obtain the amplitude in the impact parameter representation. Then, 
for given impact parameter 5o between the nuclei, the only formal difference from the hadron- 
hadron case will be the averaging over nuclear ground states, which amounts to an integration over 
transverse nucleon coordinates bf and 6j in the projectile and in the target respectively. We write 
this integration symboHcally as 

» „ A B 

/ dTAB ■■= / \{d%tTA{bt)\{d^bfTB{bf), (15) 

1=1 j=l 

with A, B being the nuclear mass numbers and with the so-called nuclear thickness function TA{h) 
being defined as the integral over the nuclear density pa(b)- 

TA{b) := J dz pAibx,by,z). (16) 

It is convenient to use the transverse distance bk between the two nucleons from the k-th nucleon- 
nucleon pair, i.e. 

bk = \bo + b^^k)-b?(k) 

where the functions 7r(fc) and r(fc) refer to the projectile and the target nucleons participating 
in the k^^ interaction (pair k). In order to simplify the notation, we define a vector b whose 
components are the overall impact parameter 60 as well as the transverse distances 61, ...,bAB of 
the nucleon pairs, 

b = {60,61, ...,6as}- 

Then the nucleus-nucleus interaction cross section can be obtained applying the cutting procedure 
to elastic scattering diagram and written in the form 

^i^fi(s) = J dHo J dTABlABis,b), (17) 

where the so-called nuclear profile function jab represents an interaction for given transverse 
coordinates of the nucleons. 

The calculation of the profile function jab as the sum over all cut diagrams of the type shown in 
fig. II does not differ from the hadron-hadron case and follows the rules formulated in the preceding 
section: 

• For a remnant carrying the light cone momentum fraction x {x~^ in case of projectile, or x~ 
in case of target), one has a factor -Fromn(a;)- 
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A 



B 




Figure 9: Example for a cut multiple scattering diagram, with cut (dashed lines) and uncut (full lines) elementary 
diagrams. This diagram can be translated directly into a formula for the inelastic cross section (see text). 

• For each cut elementary diagram (real elementary interaction = dashed vertical line) at- 
tached to two participants with light cone momentum fractions and x~ , one has a factor 
G{x~^ ,x~ , s,b). Apart from x^ and x~ , G is also a function of the total squared energy s 
and of the relative transverse distance b between the two corresponding nucleons (we use G 
as an abbreviation for G^^^ for nucleon-nucleon scattering). 

• For each uncut elementary diagram (virtual emissions = full vertical line) attached to two par- 
ticipants with Hght cone momentum fractions x~^ and x~ , one has a factor —G{x^,x~, s, 5), 
with the same G as used for the cut diagrams. 

• Finally one sums over all possible numbers of cut and uncut Pomerons and integrates over 
the light cone momentum fractions. 

So we find 




with 

n{k)—i 

-r = 1 - E ■ 

T(k)=j 

The summation indices nrik refer to the number of cut elementary diagrams and to number of 
uncut elementary diagrams, related to nucleon pair k. For each possible pair k (we have altogether 
AB pairs), we allow any number of cut and uncut diagrams. The integration variables x^ ^ refer to 
the /i'^ elementary interaction of the fc"^ pair for the cut elementary diagrams, the variables x^ 
refer to the corresponding uncut elementary diagrams. The arguments of the remnant functions 
^rcmn are the remnant light cone momentum fractions, i.e. unity minus the momentum fractions 
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of all the corresponding elementary contributions (cut and uncut ones). We also introduce the 
variables a:f'°-'and 2^^"^, defined as unity minus the momentum fractions of all the corresponding 
cut contributions, in order to integrate over the uncut ones (see below). 

The expression for ^ab sums up all possible numbers of cut Pomerons with one exception 
due to the factor (1 — ^oEmfc): one does not consider the case of all m^'s being zero, which 
corresponds to "no interaction" and therefore does not contribute to the inelastic cross section. We 
may therefore define a quantity 7ab, representing "no interaction", by taking the expression for 
7AB with (1 - ^osmj replaced by (^osmj: 



7ab(s,6) 



. AB ( h \ AB ( h ^ 

h Iab k=l {\=1 ) k=l { ^' A=l J 



A 



X 



»=1 V 7r(fc)=i / j=l \ T{k)=j 



One now may consider the sum of "interaction" and "no interaction", and one obtains easily 

7ab(s,6) +7ab(s,6) = 1. (20) 

Based on this important result, we consider jab to be the probability to have an interaction and 
correspondingly ^ab to be the probability of no interaction, for fixed energy, impact parameter 
and nuclear configuration, specified by the transverse distances bk between nucleons, and we refer 
to eq. ( po| ) as "unitarity relation". But we want to go even further and use an expansion of jab in 
order to obtain probability distributions for individual processes, which then serves as a basis for 
the calculations of exclusive quantities. 

The expansion of jab in terms of cut and uncut Pomerons as given above represents a sum of 
a large number of positive and negative terms, including all kinds of interferences, which excludes 
any probabilistic interpretation. We have therefore to perform summations of interference contri- 
butions - sum over any number of virtual elementary scatterings (uncut Pomerons) - for given 
non- interfering classes of diagrams with given numbers of real scatterings (cut Pomerons) ||2^ . Let 
us write the formulas explicitly. We have 

jAB{s,b) = V... V(l-'^oV™J / \\{\\dxljx 



f AB f rrik 

.(.,6) = E--E(i-^oE>«j / n n 

mi rriAB k=l l/J=l 

k=l I p=l ) 



k.fi k,^ 



X 

fc=l 



',s,6), (21) 



where the function $ representing the sum over virtual emissions (uncut Pomerons) is given by 
the following expression 



^AB (xP^^x'^^'^s,! 



^ AB r Ik -\ AB r h ^ 

= E---E / n n^^.A^^.^A n n n-G(^tA'*fc,A>^.&fc) 

h Iab k=l Ia=1 J k=l I A=l J 

i=\ \ 7r{k)=i I j=l \ r(k)=j I 



This summation has to be carried out, before we may use the expansion of 7a_b to obtain probability 
distributions. This is far from trivial, the necessary methods are described in ||l^. To make the 
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notation more compact, we define matrices and X , as well as a vector m, via 



X 



X+ 




m 



which leads to 



lAB{s,h) 



m 



I 



dx+dx-n^^'g'{m,x+,x-), 



lAB{s,b) 



f^iiT (0,0,0), 




n^lP{ni,X+,X-)^]l\—l[G{xl^,x-^^,s,b,) \ <i>AB{x^'°\x''^'^s,b). 



This allows to rewrite the unitarity relation eq. ( po| ) in the following form, 



/ dX+dX-n[l'g\m, X+ , X-) = 1. 



m 



This equation is of fundamental importance, because it allows us to interpret n^'^'''\m, X'^ , X^) 
as probability density of having an interaction configuration characterized by to, with the light 
cone momentum fractions of the Pomerons being given by X^ and X^ . 

7 Virtual Emissions and Markov Chain Techniques 

What did we achieve so far? We have formulated a well defined model, introduced by using the 
language of field theory, solving in this way the severe consistency problems of the most popular 
current approaches. To proceed further, one needs to solve two fundamental problems: 

• the sum ^ab over virtual emissions has to be performed, 

• tools have to be developed to deal with the multidimensional probability distribution f^^^' , 

both being very difficult tasks. Introducing new numerical techniques, we were able to solve both 
problems, as discussed in very detail in ||l0[ . 

Calculating the sum over virtual emissions is done by parameterizing the functions G 

as analytical functions and performing analytical calculations. By studying the properties of ^ab, 
we find that at very high energies the theory is no longer unitary without taking into account 
screening corrections due to triple Pomeron interactions. In this sense, we consider our work as a 
first step to construct a consistent model for high energy nuclear scattering, but there is still work 
to be done. 

Concerning the multidimensional probability distribution il^^'g (m, X^ , X^ ) , we employ meth- 
ods well known in statistical physics (Markov chain techniques) . So finally, we are able to calculate 
the probability distribution ^ab^ (to, X~), and are able to generate (in a Monte Carlo fashion) 
configurations (m, X~^ , X~) according to this probability distribution. 
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8 Summary 

What are finally the principal features of our basic results, summarized in eqs. ([l^, |2^, p^? 
Contrary to the traditional treatment (Gribov-Regge approach or parton model), all individual 
elementary contributions G depend explicitly on the light-cone momenta of the elementary in- 
teractions, with the total energy-momentum being precisely conserved. Another very important 
feature is the expHcit dependence of the screening contribution ^ab (the contribution of virtual 
emissions) on the remnant momenta. The direct consequence of properly taking into account 
energy-momentum conservation in the multiple scattering process is the vaHdity of the so-called 
AGK-cancelations in hadron-hadron and nucleus-nucleus collisions in the entire kinematical region. 

The formulas ^ |2^) allow to develop a consistent scheme to simulate high energy nucleus- 
nucleus interactions. The corresponding Monte Carlo procedure is exactly based on the cross 
section formulas so that the entire model is fully self-consistent. 
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